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ABSTRACT 


A  technique  utilizing  a  piezoelectric  disc  and  a  liquid 
filled  exponential  horn  combination  for  the  generation  and  detection 
of  ultrasound  is  described. 

The  experimental  results,  which  are  consistent  with 
theoretical  predictions,  indicate  that  these  devices  can  be 
used  to  closely  simulate  the  radiation  profile  of  point  sources 
or  receivers. 

The  data  also  suggest  that  these  transducer-horn  units 
can  be  conveniently  emnloyed  as  elements  in  the  newer  "crossed- 
array"  acoustical  holographic  systems. 
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CHAPTER  1 

INTRODUCTION 


In  recent  years  acoustics  has  enjoyed  rapidly  growing  popularity 
in  scientific  research.  The  increased  activity  in  this  area  is  largely 
due  to  the  many  new  and  imaginative  uses  that  scientists  and  engineers 
are  finding  for  ultrasonic  vibrations  [1] .  Medical  people  are  using 
echograms  as  a  method  of  forming  images  of  the  interior  anatomy. 
Materials  scientists  are  investigating  ways  of  utilizing  ultra-sound 
for  non-destructive  testing.  In  the  field  of  optical  communications 
acoustic  waves  are  used  to  amplitude  and  frequency  modulate  laser 
light.  As  a  last  example  probably  the  most  exciting  research  in  acous¬ 
tics  has  been  in  the  new  field  of  acoustic  holography. 

Because  of  the  many  uses  for  ultra-sonic  vibrations  much  attention 
has  been  directed  towards  acoustic  radiators  and  detedtors.  For 
instance  a  recent  publication  describes  the  design  of  an  ultrasonic 
generator  for  the  treatment  of  Meniere's  disease[2] .  This  thesis 
describes  a  device  which  may  be  used  as  either  a  point  source 
receiver  or  a  point  source  transmitter. 

Conventional  spherical  radiators  of  ultra-sound  in  water  are 
usually  made  of  a  piezoelectric  ceramic  which  operates  in  a  thick¬ 
ness  resonance  mode  [3] .  The  shape  of  the  radiator  which  is  that 
of  a  spherical  surface  segment  serves  two  purposes.  Firstly  the 
large  area  provides  a  method  of  coupling  a  sufficient  amount  of 
acoustic  power  into  the  water  and  secondlv  the  concave  radiating 


. 
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surface  results  in  focusing  of  the  acoustic  radiation  to  simulate 
a  point  source  radiator. 

The  device  described  in  this  thesis  consists  of  a  piezoelectric 
disc  which  radiates  into  a  horn  of  exponentially  decreasing  cross- 
section.  The  disc  provides  the  large  area  required  for  coupling 
sufficient  energy  into  the  water  while  the  horn  serves  to  focus 
the  radiation  to  closely  simulate  a  point  source.  This  device 
is  reversible  and  may  be  used  as  a  sensitive  non-directional 
point  receiver. 

In  many  respects  the  piezoelectric  disc-horn  device  de¬ 
scribed  in  this  thesis  is  superior  to  the  more  conventional 
acoustic  transducers  as  will  become  apparent  later  in  this  thesis. 
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CHAPTER  2 

DEVICE  DESCRIPTION 


2.1  Physical  details 

The  device  shown  photographically  in  figure  [1]  and  schematic¬ 
ally  in  figure  [2]  consists  of  an  acoustic  transducer  mounted  in  a 
liquid  filled  horn.  The  transducer  is  piezoelectric,  it  is  one 
centimeter  in  diameter  and  operates  in  the  thickness  mode  with  the 
resonant  frequency  being  one  megahertz.  The  horn  is  made  of  copper 
and  has  an  exponentially  decreasing  radial  cross-section.  Fabrica¬ 
tion  of  these  horns  was  achieved  by  electroforming  copper  onto  a 
stainless  steel  mandrel.  Appendix  [A]  contains  electroforming  details. 


FIGURE  1 

Photograph  of  electrof ormed  copper  horn  and  stainless  steel 


mandrel 
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Dimensions  in  centimeters 


FIGURE  [2]  pzt-horn  details 
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2.2  Qualitative  description  of  operating  principles 

The  piezoelectric  transducer  which  is  located  in  the 
large  end  of  the  copper  horn  can  perform  either  of  two  functions. 
When  the  device  is  used  as  a  transmitter  the  piezoelectric  disc 
transforms  electrical  power  into  acoustic  power  or  when  the  device 
is  used  as  a  receiver  the  transducer  delivers  an  electrical 
signal  from  the  received  acoustic  signal. 

Because  of  its  shape  the  horn  acts  as  an  acoustic  impedance 
transformer.  The  impedance  seen  by  the  crystal  is  that  of  the 
radiation  impedance  of  the  small  aperture  of  the  horn  multiplied 
by  the  ratio  of  the  transducer  area  to  the  smaller  aperture  area. 
It  is  this  impedance  transforming  property  that  enables  the  horn 
to  achieve  an  intense  acoustic  field  at  the  small  aperture  of  the 
horn. 

The  smaller  aperture  of  the  horn  illustrated  in  figure  [2] 
is,  for  the  frequency  used,  a  fraction  of  a  wavelength  in  radius. 
This  aperture  closely  simulates  either  a  point  source  of  high 
intensity  or  a  point  receiving  element  over  a  limited  field  of 
view. 

The  small  aperture  and  the  flexible  character  of  the  thin 
copper  horn  make  it  possible  to  create  an  array  of  several  point 
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receivers  very  close  together.  This  array  fabrication  is  accomp¬ 
lished  by  bending  of  the  acoustic  waves  in  a  manner  analogous  to 
the  way  in  which  light  waves  are  bent  via  fiber  optics. 

Very  recently  W.  Wells  introduced  a  new  acoustic  holographic 
system  which  he  calls  the  "multi-element  linear  crossed  array"  [4] . 
This  system  promises  to  alleviate  many  of  the  problems  found  in 
conventional  holographic  methods  [5,  6,  7], 

The  sources  and  receivers  to  be  used  in  such  a  holographic 
system  must  be  spaced  not  more  than  one  or  two  wavelengths  apart 
[8].  This  constraint  had  previously  limited  Wells’  system  to 
long  wavelength  holography  where  it  is  relatively  easy  to  make 
transducers  which  radiate  a  reasonsable  amount  of  power  from  an 
aperture  of  radius  of  less  than  one  wavelength.  The  aforementioned 
properties  of  the  piezoelectric  horn  device  meet  the  requirements 
of  the  holographic  system  proposed  by  Wells  so  that  the  crossed 
array  system  has  now  become  feasible  for  frequencies  of  one  meg¬ 
ahertz  and  perhaps  as  high  as  ten  megahertz. 


3391  lo  21/1  bsi 
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CHAPTER  3 

THEORETICAL  DISCUSSION 


3.1  Mathematical  treatment  of  horn 

The  theoretical  analysis  of  horns  which  have  exponentially 
increasing  radial  cross-sections  has  been  adequately  treated  in 
the  literature  [9,  10,  11].  Section  3.1  follows  closely  reference 
[9].  Analysis  of  horns  of  exponentially  decreasing  cross-section 
is  identical  with  the  exception  that  the  flare  constant  is  negative. 
The  cross-section  of  the  horn  is  given  by: 


S  =  S  e' 


mx 


m  = 


9  In  S 
9x 


S^  =  cross-sectional  area  at  x  =  0 
S  =  cross-sectional  area 
m  =  flare  constant 

Assuming  that  the  waves  incident  on  the  horn  are  plane  and 
that  the  phase  and  amplitude  remain  constant  on  a  plane  normal  to 
the  horn  axis  as  the  wave  propogates  down  the  horn  then  the  velocity 
potential  is  shown  in  appendix  B  to  be: 


■ 
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unx  u.x  . 

/At  1  .  „ I  2  v  lWt 

4>  =  (A  e  +  B  e  )e 


(J>  =  velocity  potential 


-.--mV 


4k^  -  =  -a  +  i3 


m 

u2  -  “  2 


4K^  -  =  -a  -  i3 


w  =  frequency  in  radians  per  second 

2tt 


k  =  wave  number  = 


X 


X  =  wavelength 


The  velocity  and  excess  pressure  are  obtained  from  the  velocity 
potential  as: 

(1)  i  =  e"ax(Ae_1(|,x  +  6)  -  Be1(fJx  +  6))eiwt 

(2)  p  =  poce"aX(Ae'1Bx  +  BeiBx)e1“t 


£  =  particle  velocity 


p  =  excess  pressure 


Po 


=  medium  density 


c  =  velocitv  of  sound  in  medium 


0  = 


-1  a 
tan  6 


Note  that  the  constants  A  and  B  are  not  the  constants  A*  and 
B'  which  appear  in  the  expression  for  .  Equations  (1)  and  (2)  are 
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written  so  that  they  may  be  easily  compared  to  plane  waves.  It 
is  noted  that  there  is  a  phase  difference  between  £  and  p  which 
does  not  occur  in  plane  waves.  Note  also  that  £  and  p  increase 
in  amplitude  exponentially  with  x. 

Knowledge  of  the  impedance  at  both  ends  of  the  horn  allows 
the  determination  of  the  constants  A  and  B.  For  a  horn  which  is 
terminated  at  x  =  0  with  a  piston  impedance  and  at  x  =  £  with 
Z  the  equations  of  force  at  the  boundaries  are: 

X/ 

(3)  at  x  =  0  Z  £  +  p  S  =  ipe^Wt 

o  1  11 

Z^  =  piston  impedance  at  the  front  of  the  horn 

,  iwt  ,  ... 

ipe  =  the  driving  force 


(4)  at  x  =  £  Z^2  ~  P2S2  =  ° 

£  =  length  of  horn 

Z  =  the  complex  impedance  which  terminates  the  horn. 

X/ 

From  equations  (1) ,  (2)  and  (4) : 


B  =  Ae 


-i2$£  /Z£e  pocS2 


7  ie  A  c 

e  +  d  cS 
£  o  2 


Now  by  using  (3)  we  have 


Z  e  +  p  cS 
A  =  { - ^ 


2De 


-i3£ 


/ 


' 


£)  ,  (I)  anoiisiipa  raoiH 
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Hence : 

B .  e-iBv  -Ae~l62;  p°c%, 


2D  = 


(7  7  i  22qq\/'  i3^/  -iBAv 

(ZoZ£  +  po  c  SlS2)(e  ‘  e  > 

+  P  c{Z  S,(ei(W  -  0)  +  e-i(^  -  6)) 
o  o  2 

+  ZtS  (e1(Bl  +  e)  +e-i(M  +  9))} 


D  =  (Z  Z  +  p^c^S  S  )i  sin[3&] 

Ox,  o  1  Z 


-f  cP  {Z  S  cos[3&  +  0]  +  Z  S_  cos[3&  -  0]} 
o  x,  1  o  2 


The  expressions  for  £  and  p  become: 

-ax 

P  ce 

p  =  •  - -  {pocS2isin[3(^  -  x)  ]  +  Z^cos[3(£  -  x)  +  0]}^e1W 

-ax 

£  ~  ~e~7T "  {Z  i  sin  [  3  (£  -  x)  ]  +  p  cS  cos[3(i  -  x)  -  0]}^ela)t 

An  equivalent  "T"  network  for  the  horn  is  developed  by 
considering  the  horn  to  have  inertia  coefficients  ,  M  and 
just  as  a  transformer. 

In  electrical  engineering,  systems  containing  four  variables 
that  are  uniquely  determined  when  two  of  the  variables  are  known 
are  called  two  port  networks.  There  are  many  ways  in  which  such 
networks  may  be  represented,  one  of  which  is  known  as  the  equivalent 
"T"  network. 


. 
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FIGURE  3 


Electrical  analog  of  horn 
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may  be  determined  by  an  open  circuit  on  the  secondary;  or 
in  the  case  of  the  horn,  =  00  so  that  ^  =  0*  L?  may  be  determined 
in  a  similar  manner.  The  mutual  impedance  =  iwM  is  the  ratio 
of  P2^2  to  ^or  closure  of  the  horn  at  S^;  i.e.,  =  °°. 

For  x  =  0  and  Z  =  00  the  expression  for  £  becomes 

X/  _L 


K,  = 


i  sin[££]i|;  e 


ioot 


1  Z  i  sin[££]  +  p  cS1  cos[££  +  9] 
o  o  1 


The  driving  point  impedance  is  then 


Z  =  Z  + 
o 


p  cS1  cos[££  +  0] 
o  1 _ 

i  sin[££] 


Subtracting  Z^  the  impedance  of  the  horn  itself  is 


iwL^  =  Z^  = 


i  p  cS, (a  -  Bcot [££] ) 
o  1 


In  a  similar  manner  the  impedance  at  the  other  end  of  the  horn 


is : 


Z2  =  icoL2  -  - 


i  p  cS„ (a  +  gcot [££] ) 
o  2 


In  finding  M  the  analysis  is  simpler  if  Z  =  0.  At  x  =  0 
with  Z^  =  °°; 


.  r  i  1  ,  i(i)t 

*  _  1  sm[££]^  e _ 

p  cS  cos[££  +  0] 
o  1 
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Again  with  Z.  =  °°  and  Z  =  0; 

Jc  o 


-a£  rQl  .  iwt 

e  cos[Q]ijj  e 
P2  S1  cos  [$£  +  0] 


The  force  at  x  =  l  is; 


eys2si  ^ elwt 


S2P2  k  S  cos[3£  +  6] 


since  e 


-ail 


S2P2 

M  =  i  M  =  -■ 

12  l 


i  P 


c  &/s7l 
o  1 


k  sin[0£] 


If  the  horn  is  driven  by  a  piston  of  negligible  impedance 
then  the  driving  point  impedance  of  a  loaded  horn  is  from  the  equiv¬ 
alent  T  network; 


Z  = 


Z1Z2  -  M12  +  hZl 

Z2  +  h 


Z  =  load  impedance 

X/ 

Z.  =  p  cS,(R  +  iX) 

R  and  X  are  frequency  dependent  quantities  which  are  dimension¬ 
less.  They  will  be  discussed  later  in  this  chapter. 

To  investigate  the  velocity  transforming  effects  of  the  horn 
we  have  from  the  equivalent  network; 
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C2)Mi2 


^2(Z2  +  V 


M 


12 


Z2  +  Z 


£ 


k2  i  cos[6]VS1S2 

^1  S2{R  sin[3£]  -  (sin[ 0] sin[3£]  +  cos[0]  cos[3£])>) 


When  3£  =irn  where  n  is  an  integer  the  velocity  transformation 
reduces  to : 


This  last  relationship  illustrates  that  the  effect  of  the 
horn  when  in  resonance  is  to  increase  the  velocity  and  pressure 
by  the  ratio  of  the  root  of  the  end  areas. 

3.2  Directivity 

In  order  to  calculate  the  pressure  distribution  at  a  distance 


. 
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away  from  the  radiating  end  of  the  horn  it  is  convenient  to  consider 
the  horn  made  up  of  many  point  sources.  A  small  element  of  area 
dS^  of  a  vibrating  source  mounted  in  an  infinite  baffle  so  that  the 
radiation  is  confined  to  a  half  space  produces  a  pressure  dp  given 
by  [12]: 


dp 


ip  ck  U.dS  el(a3t 
o  —  — _ 

2ur  ’ 


kr’) 


(see  Appendix  H) 


=  velocitv  of  the  element  dS 


FIGURE  4 

Diagram  to  obtain  directivity  function 
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In  order  to  find  the  pressure  at  a  point  in  the  half  space  it 
is  necessary  to  integrate  over  the  source  area.  The  baffle  is  assumed 
to  be  rigid  so  that  everywhere  on  the  baffle  U.dS  =  0. 


i  (u)t  -  kr '  ) 


In  the  case  of  the  horn  the  velocity  IJ  =  £  and  is  everywhere 


normal  to  S_. 

If  p  is  to  be  evaluated  at  positions  which  are  far  enough  from 
the  source  then  r’  may  be  replaced  by  r  in  the  denominator  of  the 
expression  for  p(r,0).  The  r'  which  appears  in  the  exponent  may  not 
be  replaced  by  r  since  this  term  accounts  for  the  relative  phase 
differences  due  to  radiation  from  different  positions  on  the  vibrating 
surface . 


9  9  1/9 

r’  =  {r  +  a  -  2rasin[0]  cos[i|j]} 


r*  =  r  +  ~  -  osin[0]  cos|>]  -  ^  “  ~  sin[0]  cos[^]}2  + 


In  order  that  terms  in  the  above  expression  not  contribute  to 


a 


2 


the  phase  they  must  be  <<  X.  If  r  >>  a  so  that 


r 


<<  X  then  r'  may  be 


approximated  in  the  exponent  as; 


r 


r  -  osln[0]  cos  [^] 


hence 


This  expression  is  integrated  to  obtain: 


P  = 


p  cka 
o 

2r 


Ue 


Kwt  -  kr)  2  jl(ka  sln[01) 


ka  sin[0] 


From  this  expression  we  expect  the  phase  of  the  wave  transmitted 

— - 

from  the  horn  to  be  that  of  a  spherical  wave,  i.e.,  e 

The  amplitude  is  proportional  to  r  as  in  a  sperical  wave  but 


there  is  an  additional  dependence  upon  the  azimuthal  angle.  The 


directivity  term 


2  J^(ka  sin[0]) 


is  seen  to  be  dependent  upon  the 


ka  sin[0] 

radius  of  the  radiating  aperture  and  upon  the  azimuthal  angle  0. 


3,3  Load  impedance 

The  term  Z  ,  or  the  load  impedance,  has  already  been  used  in 

X/ 

previous  sections  of  this  thesis.  This  term  is  now  discussed  in 
detail,  following  closely  reference  [13]. 


FIGURE  5 

Diagram  to  obtain  reaction  force  on  a  vibrating  piston 
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Each  element  of  a  vibrating  piston  in  an  infinite  baffle 
produces  a  pressure  wave  which  creates  a  force  on  the  rest  of  the 
piston.  A  small  element  dS  produces  a  pressure  adjacent  to  an 
element  dS1  given  by  [12], 


i  p  ck  .  ,  ,  . 

,  o  ¥T  i(wt  -  kr)  „„ 

dp  =  — -  U  e  dS 

27Tr 


Ue  =  velocity  of  the  vibrating  piston  and  is  normal  to  dS. 
For  the  horn  Ue^Wt  = 

The  total  pressure  adjacent  to  the  element  dS'  is; 


-// 


i  p  ck  .  ,  .  N 

o_  „  ei(a,t  -  kr) 

27ir 


The  total  reaction  force  on  the  piston  is  then; 


^  -  -JfP  dS- 

The  minus  sign  occurs  because  for  positive  pressures  the 
force  is  into  the  piston. 

Realizing  that  the  pressure  adjacent  to  dS’  caused  by  dS  is 
the  same  as  the  pressure  caused  by  dS'  adjacent  to  dS  makes  the 
integration  easier.  The  integration  is  carried  out  only  over  the 
inner  circle  of  radius  a  so  that  each  pair  of  elements  are  considered 
only  once.  The  result  is  then  doubled. 


dS 


r  d6  dr 


'1  -  Vi-  ; 
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The  maximum  length  in  the  direction  0  is  2ocos[0]  so  that 

7T  TT 

r  is  integrated  from  0  to  2acos[0].  The  limits  on  0  are  -  —  to  -  . 


dS'  =  adcrd^ 


The  limits  on  ip  are  0  and  2tt  while  the  limits  on  a  are  0  to  a. 


2tt  tt/2  2acos[0] 


:  =  -  U  eiwtf  ado /  d*  J  def 


dip 

0  -tt/2  '0 


-ikr  , 
e  dr 


tt/2  2acos[0] 

/  “I 

-tt/2  0 


tt/2 


-ikr  , 
e  dr 


tt _ 2i 

ik  k 


-i2ka cos[0] 
e  d0 


0 


tt/2 


l 


-i2kacos[0] 
e  d0 


tt  / 1  /01  n2  ,  1.3(2ka) 

2<1  “p-  (2k0)  +  2X4!— 


-i  ^2ko  -  (2ka)3  ..1 


The  above  integral  has  been  evaluated  by  expanding  the  expon¬ 
ential  term  and  using; 


tt/2 

|c°sn 


[x]  dx  = 


0 
tt/2 


|  cos" 
0 


[x]  dx  = 


1. 3.  .  .  (n  -  1)  ti_ 
2 . 4 .  . .  n  2 


2.4. ..(n  -  1) 

1.3. .  .  n 


n  even 


n  odd,  and  greater  than  1 


i  an  >211  fc««  0H9J  I*  L  ' 
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After  executing  the  remaining  integrations  one  obtains: 


fr  =  ~  PoC7Ta2  U  e1Wt[R(2ka)  +  iX(2ka)] 


where  R(x) 


2 

x _ 

2.4 


4 

x 

2 . 42 . 6 


+ 


6 

x _ 

2  2 

2.4  .6  .8 


•  •  • 


X(x) 


3 

x 

3.5 


+ 


5 

x 


.7 


Hence 

Zr  =  Poc7ra2[R(2ka)  4-  iX(2ka)] 

2 

For  the  horn  tt a  =  so  that 

Z  =  Z  =  p  cS  [R(2ka)  +  iX(2ka)] 
36  r  o  Z 


3.4  Design  of  horn 

The  parameters  of  the  exponential  horn  are  frequency  dependent 
so  that  as  a  starting  point  in  designing  a  horn  the  operating 
frequency  was  first  specified.  One  megahertz  was  chosen  for  the 
following  reasons: 

[1]  The  wave  length  of  one  megahertz  acoustic  waves  in 
water  is  1.48  mm.  Examination  of  the  directivity  function 


__ 


' 


-21- 


reveals  that  the  smaller  the  aperture  relative  to  the  wave 
length  the  more  widely  scattered  is  the  radiation.  A  wave 
length  of  1.48  mm.  meant  that  it  was  possible  to  machine  a 
mandrel  having  a  smaller  radius  of  the  order  of  a  wave  length. 

[2]  One  megahertz  thickness-mode  barium  titanate  piez¬ 

oelectric  discs  were  commercially  available.  The  diameter 
of  these  discs  was  one  centimeter. 

[3]  Electrical  circuitry  which  functions  at  one  mega¬ 
hertz  would  be  relatively  easy  to  build. 

The  larger  diameter  of  the  horn  was  chosen  to  be  1.5  cm.  in 
order  to  facilitate  mounting  of  the  1  cm.  crystal. 

In  designing  the  smaller  end  of  the  horn  three  factors  had 
to  be  considered. 

[1]  Directivity  function 

[2]  Radiation  impedance;  R(x)  +  iX(x) 

[3]  Ease  of  machining  a  mandrel  having  a  sufficiently 
small  radius  at  the  smaller  end. 

By  choosing  the  smaller  diameter  of  the  horn  to  be  1.5  mm. 
(approximatelv  one  wave  length  at  one  megahertz)  the  directivity 
criterion  would  be  satisfied,,  At  an  azimuthal  angle  of  15  degrees 
the  amplitude  is  96  percent  of  the  on  axis  amplitude.  The  values 
for  R(x)  and  X(x)  are  1  and  .,2  respectively  so  that  an  appreciable 
amount  of  the  power  that  is  incident  on  the  apperture  is  transmitted. 


- 


' 
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If  the  radius  of  the  aperture  becomes  much  less  than  X  then  the 
terms  R(x)  and  X(x)  approach  0.  For  small  apertures  (r  <<  X) 
high  acoustic  fields  are  necessary  to  produce  adequate  radiated 
power.  Since  high  fields  cause  degassing  and  cavitation  they 
are  undesirable.  Even  for  the  radius  used  occasionally  bubbles 
appeared  in  the  narrow  neck  of  the  horn  indicating  degassing.  An 
air  bubble  in  the  end  of  the  horn  greatly  deteriorates  the  horn 
performance  as  is  illustrated  in  chapter  4. 

Only  the  flare  constant  "m"  remains  to  be  specified.  Factors 
affected  by  this  parameter  are  the  cut-off  frequency  and  the  phase 
angle  between  pressure  and  velocity.  The  terms  u^  and  u^  which 
appear  in  the  expression  for  the  velocity  potential  become  entirely 
real  for  certain  values  of  m. 

u^  =  -a  +  iB 
u2  =  -a  -  13 


Hence  the  waves  become  evanescent  for  m  =  2k  =  - — 


The  cut-off  frequency  is  defined  as  the  frequency  for  which 
m  =  2k. 


A  value  of  -.848  cm 


-1 


was  chosen  for  m  so  that  the  cut-off 


frequency  is  specified  as  10  k  hertz. 


For  a  tuned  horn  (8£  -  nff ; 
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n  is  integer)  the  phase  angle  between  pressure  and  velocity  is 
given  by  0. 

.  -1  ot 

9  =  tan  j 

0  <  1  degree 

The  small  value  of  0  results  in  a  power  factor  very  close  to 

one . 
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CHAPTER  4 

EXPERIMENTAL 


4.1  Tank  construction 

Mounting  details  for  the  transmitting  horn  and  detector 
are  presented  in  figure  [6] .  The  water  tank  which  measured 
2x1x1  meters  was  constructed  of  3/4  inch  plywood.  The  size 
of  the  tank  was  larger  than  the  scanning  apparatus  required.  The 
large  tank  size  resulted  in  a  path  difference  between  signals 
reflected  from  the  walls  of  the  tank  and  the  signals  of  interest. 
Because  the  experiments  were  run  on  a  pulsed  basis  the  path 
difference  resulted  in  isolation  of  reflected  signals  from  signals 
of  interest.  The  horn  was  mounted  in  the  center  of  one  of  the  ends 
of  the  tank  so  that  the  small  aperture  would  be  in  an  infinite 
baffle  and  so  that  the  back  of  the  piezoelectric  disc  would  not 
be  mass  loaded  by  water. 

In  order  to  reduce  any  electrical  pick  up  through  the  water 
the  side  of  the  crystal  facing  the  water  was  grounded. 

Another  piezoelectric  disc  was  mounted  so  that  it  would 
rotate  about  the  aperture  of  the  horn.  A  linear  potentiometer 
was  mechanically  fastened  to  the  rotating  arm  so  that  it  was  pos¬ 
sible  to  obtain  an  electrical  signal  proportional  to  the  azimuthal 
angle  e  which  indicates  the  position  of  the  plane  mounted  disc 


relative  to  the  horn  aperture. 


in b 5  -»ri : 
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Rota  ting  Arm 


FIGURE  6 
Tank  geometry 

The  rotating  arm  is  such  that  the  detector  rotates  about  the 
horn  aperture  in  a  plane  parallel  to  the  bottom  of  the  tank. 
The  horn  axis  lies  in  the  plane  of  rotation  of  the  disc. 


* 

. 
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4.2  Radiation  pattern  for  single  horn 

4.2.1  Amplitude  versus  azimuthal  angle 

In  this  section  experimental  measurement  of  the  directivity 
function  of  the  horn  device  is  described.  The  effects  of  various 
mounting  factors  on  directivity  are  illustrated  by  means  of  plots 
of  amplitude  versus  azimuthal  angle. 

The  horn  was  used  as  a  transmitter  while  the  plane  mounted 
crystal  functioned  as  a  receiver.  Electrical  driving  signals 
were  one  megacycle  sine  waves  of  forty  volt  amplitude  while  the 
received  electrical  signals  were  on  the  order  of  a  millivolt. 

This  meant  that  electrical  isolation  of  100  db.  or  greater  was 
required  in  order  to  obtain  meaningful  measurements.  Pulsing 
the  acoustic  transmitter  provided  adequate  electrical  isolation 
of  the  receiver.  Wave  trains  of  about  .1  milliseconds  (100  cycles) 
spaced  .9  milliseconds  apart  were  used  so  that  the  acoustic  wave 
train  reached  the  receiver  in  the  absence  of  stray  electrical 
driving  signals  (the  ratio  of  stray  electrical  signals  from  the 
transmitter  to  noise  from  receiver  was  <  1  ) .  Another  positive 
effect  of  making  measurements  on  a  pulsed  basis  was  that  acoustic 
wall  reflections  and  standing  waves  in  the  tank  were  minimized. 

Received  signals  were  amplified  by  a  low  noise  amplifier 
which  was  specially  designed  to  work  at  megacycle  frequencies. 
Further  processing  involved  a  circuit  which  measured  the  am¬ 
plitude  of  the  wave  trains  and  held  this  value  until  the  next 


fjj  'ini  'B  1 1  1  cA 
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FIGURE  7 

Schematic  of  experimental  set  up  for  measuring  directivity 


function 
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FIGURE  8 


Photograph  of  experimental  equipment 
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wave  train  was  presented  to  it. 

A  circuit  diagram  of  the  gated  oscillator  amplifier  which 
provided  the  driving  electrical  signals  appears  in  appendix  C. 

The  low  noise  preamplifier  and  amplitude  detector  schematics  also 
appear  in  appendix  C. 

In  section  3.2  the  theoretical  directivity  function  for 

a  vibrating  piston  in  an  infinite  baffle  was  shown  to  be 

2  (k  a  sin[©])  .  Experimental  mounting  conditions  which  produced 
k  a  sin[e] 

results  most  closely  approximating  the  theoretical  predictions  are 
listed  below. 

[1]  The  horn  aperture  was  mounted  tightly  in  the 
baffle . 

[2]  The  wall  thickness  of  the  horn  was  as  thin  as 
possible. 

[3]  The  piezoelectric  disc  was  mounted  perpendicular 
to  the  axis  of  the  horn. 

Profiles  obtained  under  the  conditions  mentioned  above  were 
more  directional  than  the  theory  predicted.  This  is  illustrated 
in  figure  [10]  .  Such  an  effect  is  equivalent  to  a  larger  effective 
radiating  apperture.  In  order  to  develop  a  theory  which  would 
predict  accurately  the  directivity  function  of  the  horn  a  computer 
program  was  set  up  which  would  assign  a  fraction  of  the  power 


' 

. 
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radiated  to  the  copper  walls  of  the  horn  and  the  remainder  to 
the  water  filled  aperture. 


FIGURE  9 

Radiating  aperture 

The  amplitude  of  vibration  within  the  water  filled  area  is 
normalized  to  1  while  the  amplitude  of  vibration  in  the  copper 
area  is  assumed  to  be  some  fraction  "f"  of  that  of  the  center 
area.  Assuming  linear  superposition  the  amplitude  at  a  large  dis¬ 
tance  frcm  the  horn  is  proportional  to  J1  (k  r  sin[e]) 

1  w 

k  r  sin[e] 
w 

for  the  center  area  and  is  proportional  to 
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f  x 


J-j  (k  rcu  sin[0])  J1  (k  r_sin[0] 


w 


k  r  sin[0] 
cu 


k  r  sin[0] 
w 


for  the  ring  shaped  copper  area.  The  total  directivity  function 


is  proportional  to: 


J. (k  r  sin[ 0] ) 
1  _1 _ cu 

k  r  sin[0l 
cu 


'+  (1  -  f) 


J  (k  r  sin [0] 
1  w 

k  r  sin[0] 
w 


The  APL  program  presented  in  Appendix  D  was  written  to 

evaluate  the  above  equation.  Different  values  of  f  were  used 

and  then  the  calculated  results  were  compared  to  the  experimental 

results.  Closest  correlation  occurred  for  f  =  1.  Hence  for  cal- 

2  J^(k  a  sin[0]) 

culations  the  directivity  function  - ; - ; — rr-z - where  a  is 

k  a  sin[0] 

the  outer  diameter  of  the  horn  should  be  used. 

The  preceding  results  strongly  suggest  that  the  metal 
horn  itself  vibrates.  This  result  is  not  surprising  when  one 
considers  the  lack  of  rigidity  of  the  thin  electroplated  structure. 

At  this  point  it  should  be  noted  that  the  directivity  function 
above  assumes  that  the  velocity  of  the  water  aperture  is  in  phase 
with  the  velocity  of  the  copper.  The  agreement  between  experiment 
and  calculations  using  f  =  1  suggests  that  acoustic  coupling  to 
the  horn  shell  is  not  a  consequence  of  the  0  ring.  If  the  0  ring 
were  the  driving  mechanism  then  the  vibrating  copper  at  the  aperture 
would  be  out  of  phase  with  the  vibrations  in  the  aperture  since  the 
velocity  of  sound  in  copper  is  much  different  than  in  water. 


. '  . 


. 


' 


-32- 


By  the  theorem  of  reciprocity  [14]  the  horn  should  display 
the  same  directional  characteristics  when  used  as  a  receiver  or 
as  a  transmitter.  Experimentally  this  was  found  to  be  the  case. 

The  plane  mounted  disc  was  used  as  the  transmitter  and  rotated 
about  the  horn  which  was  used  as  a  receiver.  In  this  manner 
plane  waves  were  presented  to  the  horn  aperture  over  a  range  of 
azimuthal  angle  of  -  30  degrees.  Normalized  plots  of  amplitude 
versus  azimuthal  angle  were  the  same  for  both  transmitter  and 
receiver.  This  result  again  indicates  that  the  larger  effective 
aperture  is  due  to  vibration  of  the  copper  walls.  Another  result 
worthy  of  note  is  that  the  directionalitv  does  not  increase  with 
stronger  acoustic  fields  but  is  constant  over  a  wide  range  of 
intensities . 

Various  different  mounting  conditions  that  were  tried  are 
illustrated  graphically  in  figures  [12]  thru  [16] .  Since  the  metal 
of  the  horn  was  thought  to  be  vibrating  it  was  possible  that  by 
contouring  the  end  of  the  horn  less  directivity  could  be  obtained. 
However  the  best  results  were  obtained  by  simplv  making  the  walls 


of  the  horn  as  thin  as  possible. 
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DIRECTIVITY 


FIGURE  10 


Comparison  of  theoretical  and  calculated  directivity  for  horn 
aperture  of  inner  radius  .1  cm.  and  outer  radius  of  .15  cm. 
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DIRECTIVITY 


FIGURE  11 

Comparison  of  experimental  data  with  modified  theory  which 
takes  into  account  vibrating  copper.  (For  modified  theory 
see  pages  2  9  -  31  .) 
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DIRECTIVITY 


Azimuthal  Angle  (degrees) 


FIGURE  12 


Directivity  for  thick  walled  horn 
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EXPERIMENTAL 

DIRECTIVITY 


Azimuthal  Angle  (degrees) 


FIGURE  13 


Directivity  for  concave  horn  end 
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EXPERIMENTAL 

DIRECTIVITY 


Directivity  for  horn  with  concave  end  protruding  from  baffle 
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EXPERIMENT  AL 
DIRECTIVITY 


Directivity  of  horn  with  convex  end 
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EXPERIMENTAL 

DIRECTIVITY 


FIGURE  16 


Convex  horn  end,  baffle  aperture  filled  with  wood  putty 
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EXPERIMENTAL 

DIRECTIVITY 


FIGURE  17 

Comparison  of  directivity  plots  for  horn  used  as  a 


transmitter  and  as  a  receiver 
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EXPERIMENT  AL 
DIRECTIVITY 


Directivity  for  horn  with  air  bubble  in  aperture. 

Due  to  normalization  the  plot  does  not  show  the  greatly  reduced 
signal  level.  Signal  levels  were  so  small  that  electrical  noise 
prevented  obtaining  a  reasonable  plot. 
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FIGURE  19 


Photograph  of  various  radiating  apertures 
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FIGURE  2  0 

Oscilloscope  trace  illustrating  electrical  isolation  of 
receiver  and  transmitter.  The  upper  trace  is  the  electrical 
signal  applied  to  the  transmitter  while  the  lower  trace  is 
the  received  signal.  Time  separation  between  the  two  signals 
is  due  to  the  distance  between  the  transmitter  and  the  receiver 
and  the  velocity  of  sound.  Echoes  from  the  walls  of  the  tank 
are  apparent  on  the  lower  trace. 
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4.2.2.  Phase  versus  azimuthal  angle 

For  large  distances  from  the  radiating  aperture  the  pressure 
is  given  by: 


P 


i  p  a  ck  ,  . 

- 2 -  u  "  kr) 

2r 


2  J^(k  a  sin[0]) 
k  a  sin[0] 


As  was  mentioned  in  chapter  3  this  equation  illustrates 
that  the  phase  of  the  radiation  is  expected  to  be  that  of  a 
spherical  wave,  (e  ) .  In  order  to  test  this  two  experiments 
were  carried  out.  The  first  arrangement  is  illustrated  in  figure 
[21]  and  is  fairly  straightforward.  The  received  acoustic  wave 
was  displayed  on  one  channel  of  the  CRT  while  the  electrical  wave 
used  to  drive  the  transducer  was  displayed  on  the  other  channel  to 
provide  a  reference.  The  receiving  element  was  rotated  about  the 
transmitting  horn  at  a  constant  radius  while  the  waves  on  the  CRT 
were  monitored.  No  phase  shift  was  observed  indicating  that  the 
phase  dependence  of  the  radiation  is  that  of  a  spherical  wave. 

However,  the  mechanical  scanning  structure  was  not  stable 
so  that  another  experiment  involving  two  transmitters  was  set  up 
in  order  to  further  examine  the  phase  dependence  of  the  radiation 
from  the  horn.  The  basic  geometry  of  the  experiment  which  is  shown 
in  figure  [22]  is  not  critically  dependent  upon  r.  Unfortunately 
the  results  from  this  experiment  are  also  dubious  for  reasons 


■ 
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FIGURE  21 

Experimental  set  up  for  first  experiment  to  measure 


phase  dependence 
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FIGURE  2  2 

Geometry  for  the  second  experiment  to  measure 
phase  dependence 
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which  will  be  explained  latter.  This  experiment  is  presented  to 
illustrate  some  of  the  problems  in  making  a  phase  measurement  on 
the  horns. 

This  experiment  is  similar  to  Young's  experiment  in  that  an 
interference  pattern  from  two  sources  is  recorded.  However,  the 
recording  of  the  interference  pattern  is  made  over  an  arc  instead 
of  on  a  planer  surface.  In  addition  the  recording  is  made  over 
a  much  wider  angle  so  that  approximating  the  waves  to  be  plane 
is  not  valid  as  in  Young's  experiment. 

Referring  to  figure  [22]  the  values  of  and  are  given 

by: 

(r ,  0)  =  [  (r  cos[0])2  +  (r  sin[0]  +  d)2]1//2 

L2  =  L2(r,0)  =  [ (r  cos[0] )2  +  (r  sin[e]  -  d)2]1/2 


For  and  L >>  a  then  the  amplitude  profiles  for  horns  one 
and  two  are  given  by: 

Af  (k  a  sin[0^]  ) 

1  L^k  a  s  in  [  0  ^  ] ) 

A?  J,  (k  a  sin [0 „ ]  ) 

D  =  - 

L>2  (k  a  sin[©2] ) 

Let  the  phase  dependence  be  given  by: 

e-i(kLl  +  f.ce^)  .  ei+x 
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l2  +  f202) 


2 


The  time  dependence,  e1Wt,  is  not  included  here  as  it  is  not 
necessary  for  the  analysis.  The  above  expressions  refer  to  horns 
one  and  two  respectively.  Note  that  the  phase  dependence  is  that 
of  a  spherical  wave  if  the  functions  f^(0  )  and  f 2 ( ® 2 )  are  zero* 

When  and  are  sufficiently  large  so  that  the  expressions 
for  and  are  valid  then  the  phasor  (Ae1^)  of  the  radiation 
formed  by  the  interference  of  the  radiation  from  the  two  horns  may 
be  found  by  referring  to  figure  [23] 


FIGURE  23 

Phasor  diagram  used  to  calculate  the  interference  pattern  from 


two  horns. 


. 
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Ae1'*’  =  D1e1<t’l  + 

A  =  |Ae1<P|  =  {(D1  +  D2cos[tj)2  -  1  ] )  +  (D2sin[(}>2  -  tj^])  } 

It  is  this  function  that  the  experiment  attempts  to  measure, 
and  L2  and  r  have  approximately  the  same  direction  0  if 
the  position  of  the  detector  in  space  is  such  that  r  >>  d.  This 
being  the  case  then  0^  E  02  H  6 •  and  L2  may  be  replaced  by  r 

in  the  denominator  of  the  expressions  for  and  D2  if  r  is 
sufficiently  large.  Hence: 

Af  J^(k  a  sin[0^])  A^  J^(k  a  sin[0]) 

D1  L^(k  a  sin[0^])  r (k  a  sin[0]) 

A2  (k  a  sin [ @2 ] )  J^(k  a  sin[0]) 

D  L2(k  a  sin[02])  r(kasin[0]) 

If  the  radiation  from  the  horns  has  spherical  phase  dependence 
then  f^(0^)  anc*  ^2^2^  =  0  so  that  the  expresssions  for  and  4> 2 
become : 

2ttL 

*i  =  kLi  ‘  — 


2nL, 

4  2 =  kL2  “  ir 


Hence : 


A  = 


(k  a  sin [0 ] ) 
r  (k  a  sin[ 0 ] ) 


A1  +  A2cos(  (L2  -  Lp-^] 


+ 


A2sir.  (  (L2  -  Lx)^. 


1/2 


L 


' 
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Experimental  conditions  correspond  to  holding  r  =  R  =*26.6  cm. 
=  constant  and  varying  0  from  -30  to  +30  degrees.  The  function 
A  was  evaluated  for  these  conditions  by  means  of  the  computer 
program  in  appendix  (E) .  A  comparison  of  the  experimental  and 
calculated  results  appears  in  figure  [24] .  The  difference  in  the 
amplitude  envelopes  may  be  attributed  to  the  directivity  functions 
varying  slightly  from  their  expected  values.  A  more  serious 
discrepancy  between  the  experimental  and  theoretical  curve  is 
that  for  increaseing  0  the  locations  of  the  minima  and  maxima  are 
not  the  same  for  the  two  curves. 

These  discrepancies  in  position  of  maxima  and  minima  can 
not  be  attributed  to  variations  in  R  from  its  nominal  value  of 
26.6  cm.  As  was  previously  mentioned  this  experiment  is  not 
critically  dependent  upon  R  and  this  is  indeed  the  case.  By 
means  of  the  computer  the  difference  described  below  was  calculated. 

A(r=26 .6,0)  _  A(r=30 , 0 ) 

A(r=26.6,0)  A(r=30,0) 

-4 

The  greatest  difference  was  less  than  2  x  10 

Since  the  horns  are  very  similar  in  construction  it  seems 
reasonable  to  assume  that  the  phase  dependence  for  both  horns  is 
the  same.  In  other  words  it  would  be  reasonable  to  assume  that 
f 1 (0 )  =  f  ^  ( 0 ) .  Providing  that  f  x  (0 )  does  not  vary  rapidly  with 


. 
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INTERFERENCE  PATTERN 


Interference  pattern  obtained  from  two  radiating  horns,  compared 
to  the  interference  pattern  from  two  sources  having  spherical 
wave  phase  dependence.  The  radius  of  the  arc  the  detector 


travelled  was  R  =  26.6  cm. 


. 


-52- 


0  then  f^(0^)  =  -^2^2^  =  because  of  the  condition  that 

0  -  0  -  0^.  However  this  variation  from  spherical  phase  dependence 
does  not  explain  the  difference  in  locations  between  the  maxima  of 
the  two  curves.  Examination  of  the  expression  for  A  shows  that  A 
is  a  function  only  of  the  phase  difference  between  the  two  waves. 
Hence  if  f^(0^)  =  f  2  ( ©  2 )  ’  t^ien  bhese  functions  do  not  affect  A. 

One  possible  explanation  is  that  due  to  small  differences  in 
the  horns  the  functions  f  (0 )  and  are  different,  i.e., 

f^(0)  -  f2(0)  ^  0.  For  the  sake  of  illustration  let: 

f  2  (0)  -  f]_(0)  =  f  (9)  =  -.85  sin  [  0  ] 

Under  these  conditions  the  normalized  function  for  A  is 
written  as: 


A  plot  of  this  curve  appears  in  figure  [25]  and  is  compared 
to  the  experimental  curve.  Better  agreement  between  the  two  curves 
in  figure  [25]  than  in  figure  [24]  is  noted. 

Another  possible  source  of  error  is  that  the  scanning 
arc  does  not  lie  in  the  plane  defined  by  the  horn  axis  and  the 
line  joining  the  two  apertures.  Assuming  the  scanning  arc  to 
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INTERFERENCE  PATTERN 


Interference  pattern  obtained  from  two  radiating  horns  compared 
to  two  sources  having  the  phase  dependence  described  on  page  [52] 
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x 


FIGURE  26 

The  geometry  here  illustrates  the  conditions  for  the  "calculated" 
curve  anpearing  in  figure  [27],  The  plane  containing  the 
arc  of  the  detector  is  displaced  by  an  amount  X  =  5  cm.  from 
the  line  joining  the  two  sources. 


' 
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INTERFERENCE  PATTERN 


Interference  pattern  obtained  from  two  radiating  horns  compared 
to  the  expected  pattern  under  conditions  illustrated  in  figure  [26]. 
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be  defined  as  in  figure  [26]  then  the  associated  interference 
pattern  would  be  as  in  figure  [27],  The  phase  dependence  has  been 

assumed  to  be  spherical  f(0)  =  f(0)=O.  L  and  L_  in  this  case  are 

12  12 

given  by: 


[ (R  cos [ 6 ] ) 2 
[ (R  cos [ 0  ] )  2 


+  (R  sin [ 0 ]  +  d)2 
+  (R  sin[0]  -  d)2 


+  x2)] 
+  x2)] 


1/2 

1/2 


The  function  A  was  evaluated  using  these  expressions  for 
and  .  The  values  for  X  and  R  were  X  =  5  cm  and  R  =  26.6  cm. 

A  comparison  of  the  experimental  curve  with  the  function  A  is 
given  in  figure  [27].  If  compared  to  figure  [24]  an  improvement 
in  correlation  between  the  experimental  and  theoretical  curve  is 
noted . 

Choosing  X  to  be  5  cm  is  a  bit  extreme  since  the  mechanical 
scanning  is  more  accurate  than  that.  With  this  in  mind  it  is 
reasonable  to  assume  that  the  phase  dependence  of  the  two  horns 
is  not  that  of  a  spherical  wave,  nor  is  the  phase  dependence  of 
the  two  horns  the  same  since  it  has  been  shown  that  the  discrepancies 
which  appear  in  figure  [24]  can  be  attributed  to  f^CQ)  -  f (0 )  ^  0. 

If  the  radiation  characteristics  of  one  of  the  sources  were 
completely  known  (amplitude  and  phase)  the  radiation  from  the 
other  could  be  uniquely  determined  from  the  interference  pattern. 


* 

_ 
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In  other  words,  if  one  of  the  horns  could  be  made  to  interfere 
with  a  plane  wave  emitted  from  a  large  flat  crystal  then  the 
phase  characteristics  could  be  determined.  Another  alternative 
to  obtaining  accurate  phase  measurements  is  with  the  use  of  a 
scanning  apparatus  with  sufficient  stability  so  as  to  render  the 
results  of  the  first  experiment  accurate. 


4.3  Horn  and  plane  mounted  transducer  profile 

A  comparison  of  the  radiation  profile  given  by  the  horn  device 
and  that  produced  by  a  standard  plane  mounted  piezoelectric  crystal 
is  shown  in  figure  [28].  The  improvement  in  performance  character¬ 
istic  is  clearly  illustrated  in  this  diagram. 
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DIRECTIVITY 


Azimuthal  Angle  (degrees) 


FIGURE  28 

A  comparison  of  directivity  functions  for  a  plane  mounted 
crystal  and  a  horn 


' 
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4.4  Amplitude  versus  frequency 

In  chapter  3  an  expression  for  the  driving  point  impedance 
was  given  as : 


Z  = 


Z1Z2  -  M12  +  hh 


zn  +  z 


£ 


If  the  values  for  Z^,  Z^,  Z^  and  are  substituted  then: 

Pq2c^S1S2  +  i|  cS^I (m/2)  -  gcot [ 3&] }{pQcS2 [R(x)  +  iX(x)]} 
ip  cS 

- r -  ((m/2)  +  Bcot[3A])  +  p  cS  [R(x)  +  iX(x)] 

K  O  l 

Note  that  this  is  the  impedance  presented  to  the  crystal. 

Since  the  driving  point  impedance  is  a  function  of  frequency 
we  expect  the  transmitted  intensity  should  also  be  a  function  of 
frequency.  Such  an  effect  should  be  manifest  as  a  fine  structure 
superimposed  on  the  normal  "Q"  curve  of  the  crystal. 

An  investigation  of  this  effect  was  carried  out  by  means  of 
the  experiment  schematically  illustrated  in  figure  [29].  Plots  of 
transmitted  amplitude  versus  frequency  for  the  crystal-horn  device 
and  for  a  plane  crystal  were  obtained.  These  curves  appear  in 
figure  [ 30 ] . 

Also  appearing  in  figure  [30]  is  a  calculated  curve  for 
amplitude  versus  frequency.  This  curve  was  obtained  by  evaluating 


-■ 
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the  real  part  of  the  function  for  Z,  ie.  the  real  part  of  the 
impedance  presented  to  the  crystal.  This  result  was  then  multiplied 
by  the  envelope  of  the  experimental  curve  which  appears  in  figure  [30] . 
The  computer  program  used  to  carry  out  the  calculation  appears  in 
appendix  (F) . 

The  spacing  between  the  amplitude  peaks  is  determined  by  the 
term  3£.  An  increment  of  tt  in  (3£  corresponds  to  the  spacing  between 
two  adjacent  amplitude  peaks.  These  results  indicate  that  the  horn 
may  be  ’tuned'  simply  by  translation  of  the  transducer  slightly 
along  the  horn  axis.  The  ability  to  tune  each  structure  is  a  highly 
desirable  feature  from  a  system  point  of  view  as  it  means  that  only 
one  master  driving  oscillator  is  necessary  for  many  horns. 
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WT 


FIGURE  2  9 

Schematic  of  experiment  to  obtain  frequency  versus  amplitude 


plot 


-63- 


4.5  Characteristics  of  horn  in  an  array 

The  horns  may  be  bent  to  facilitate  close  mounting.  This 
is  shown  in  figure  [31].  The  next  figure  illustrates  the  effect 
on  the  radiation  pattern  due  to  close  mounting.  An  increase  in 
directivity  over  that  of  a  single  horn  is  observed  although  for 
fields  of  view  of  less  than  15  degrees  the  horns  can  be  used  to 
closely  simulate  an  array  of  point  sources  or  point  receivers. 
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FT CURE  31 

(’lose  mounting  details  to  facilitate  simulation  of 


Doint  source  array 
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DIRECTIVITY 


Azimuthal  Angle  (degrees) 


FIGURE  32 


Directivity  plot  for  horn  mounted  in  an  array 
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CHAPTER  5 

CONCLUSIONS 


The  experimental  results  presented  in  this  thesis  indicate 
that  the  PZT-horn  devices  perform  closely  to  the  expected  theoret¬ 
ical  predictions.  They  may  be  used  effectively  as  a  point  source 
of  high  intensity  or  as  a  receiver  with  a  wide  field  of  view. 

The  flexible  nature  of  the  horn  facilitates  close  mounting 
of  apertures  and  hence  the  fabrication  of  an  array  of  point 
sources  or  point  receivers. 
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CHAPTER  6 

SUGGESTIONS  FOR  FURTHER  WORK 


The  very  positive  experimental  results  obtained  from  the 
horn  described  in  this  thesis  warrant  further  investigation  into 
inverse  horns . 

Different  flare  constants  and  aperture  sizes  should  be  tried. 

A  conical  horn  would  be  much  easier  to  fabricate  than  an  exponential 
horn  and  should  be  investigated. 

Materials  other  than  copper  could  be  used.  Perhaps  a  more 
rigid  material  would  stop  the  horn  aperture  from  vibrating  and  a 
wider  field  of  view  would  be  obtained.  If  horns  were  fabricated 
by  glass  blowing  then  very  small  apertures  could  be  obtained  and 
the  horn  would  be  effective  as  a  point  source  at  higher  frequencies. 

This  thesis  has  shown  that  it  is  possible  to  use  the  horn  in 
an  array  so  that  is  is  feasible  to  realize  the  linear  crossed-array 
holographic  system  proposed  by  Wells.  Such  a  system  is  presently 
under  construction  at  the  University  of  Alberta. 
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APPENDIX  A 

ELECTROFORMING  DETAILS 


A.l  Mandrel  preparation 

For  cleaning  and  easy  removal  of  electroformed  horns  from 

the  stainless  steel  mandrel  the  mandrel  was  treated  cathodically 

in  a  6  oz,  /  gal.  solution  of  trisodium  phosphate  in  water. 

The  temperature  of  the  solution  was  between  140  and  180  degrees 

2 

Farenhiet.  A  current  of  50  amp  /ft  for  one  minute  was  used 
and  then  the  mandrel  was  rinsed  in  water  at  100  degrees  F. 

A, 2  Electroforming  solution 

The  electroforming  solution  consisted  of  30  oz .  /  gal.  of 

copper  sulfate  and  6  oz.  /  gal.  of  sulfuric  acid.  A  current 

2 

density  of  100  amps  /  ft  or  about  .5  amps  per  mandrel  was  used. 
The  solution  was  maintained  at  a  temperature  of  80  degrees  F  and 
was  constantly  circulated  about  the  mandrel. 


. 
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APPENDIX  B 
VELOCITY  POTENTIAL 


B.l  Preliminaries 

This  appendix  follows  closely  reference  [15] . 

For  small  increments  of  density  a  term  called  condensation 
(s)  is  defined  as: 


p  =  p  +  6p  =  p  (1  +  s)  (b.l) 

o  o 

since  pv  =  p  v 
o  o 

we  have  s  -  -  —  (b.2) 

v 

o 

v  =  specific  volume 
p  =  density 

v  =  mean  specific  volume 
o 

p  =  mean  density 
o 

For  an  adiabatic  fluid  which  behaves  as  described  by: 

p.vy  =  P  VY 
11  o  o 

P^  =  pressure  at 

P  =  initial  pressure 
o 

V  =  volume  at  P^ 

V  =  initial  volume 
o 


- 
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the  following  relationship  holds: 

6P  y6v 

_  =  -  —  =  YS 

o  o 

Defining  k' ,  coefficient  of  cubic  elasticity,  as 

(SP 

k’  =  YPo  =  ™  (b . 3) 

o  s 

then  P1=P  +k's=P  +6P 
1  o  o 

B.2  Equation  of  continuity 

Next  the  equation  of  continuity  is  derived. 


Considering  flow  in  the  x  direction  through  the  volume 
element  the  expression  for  the  net  flux  is: 


6 1 


pi  -  ( pi  +  3g"^dx) 


.  ,  3(pO  ^ 

A  =  -  Adx  — - -  d  t 

dX 


where  A  =  dy  dz 


Kjl  ■■  H  *  K 
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The  increase  in  the  amount  of  fluid  contained  in  the  element 
is  given  by: 


6t  A 


9p 

dx  37 


hence  the  equation: 


3p  ,  3(p|) 

9t  9x 


(b .  4) 


3  (  A  ^  ) 

using  p  =  p  (1  +  s)  and  realizing  that  —  << 

O  o  X 


it 

9x 


then : 


9s  H 

9 1  9x 


=  0 


(b.5) 


The  above  is  the  equation  of  continuity. 

•  • 

The  validity  of  assuming  — — —  <<  —  is  verified  later 
by  equations  (b.25)  and  (b.26). 

From  (b.5)  we  obtain: 


2  2* 

38  +^4=o 


9t9x 


9x 


(b .  6) 


B.3  Force  equation 

Equating  the  net  force  on  the  element  to  the  rate  of  change 
of  momentum  of  the  element  in  the  x  direction  one  obtains: 


dx  dy  dz  +  p£)  =  0 


(b .  7) 


Y^ibilav  sift 


. 
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Here  p  is  the;  excess  pressure ;  (SI’  =  p  =  k  s 


(b .  8) 


Substituting  for  p  the  equation  for  force  becomes 


A  •  • 

k*  ~  +  pt,  =  0 
9x 


(b.9) 


The  next  step  is  differentiation  with  respect  to  time  and 

•  •  • 

•\  r  r  ^  c 

the  use  of  the  approximation  that  - —  >>  - -  . 

9 1  9x 

This  approximation  is  considered  in  reference  [16]  and  is 

valid  for  waves  in  which  —  <<  1.  Referring  to  equation  (b.24) 
t  -4 

we  see  that  —  <  10  <<  1 

c 


2  2* 

k.  9  s  U  _ 

k  9x8t^\  2  "  ° 

9 1 

By  comparison  of  this  equation  to  (b.6)  the  result  is  the 
wave  equation. 


2  •  2  • 

9  £  k’  91 


(b .  10) 


9 1‘ 


9x 


u  2  k’ 

Hence  c  =  — 


(b .  11) 


where  c  is  the  velocity  of  propagation, 


B.4  Velocity  potential 

Previously  the  expression  for  forces  in  the  x  direction  was 
given  as: 
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Using  p  =  k’s  and  c  =  k'/p  expressions  for  the  three  directions 
are  obtained  as: 


Z  = 
h  = 


2  3s 
c  — 
dX 


2  3s 


Z  =  - 


2  3s 


(b . 12) 

(b .  13) 

(b . 14) 


To  simplify  the  theory  a  single  function  called  the  velocity 
potential  (cf>)  which  contains  all  three  displacement  quantities 
C,  n  and  Z  is  sought.  For  a  fluid  in  which  there  is  no  circulation 
a  function  (f)  exists  such  that: 


c 


3£_ 

3x 


n  =  - 


3(f) 

3y 


4  = 


_  ii 

3z 


(b . 15) 
(b . 16) 
(b .  17 ) 


or  <f>  =  -  j  (£dx  +  ndy  +  £dz) 

Integration  of  the  equations  (b.12),  (b.13)  and  (b.14)  yields 


Z  =  - 


n  =  -  c 


a 

3x  J 


sdt  +  Z 


2  1 

3y 


5  =  -  c 


2  3 


3z 


sdt  +  n 


o 


Sdt  +  £ 


where  L  .  n  and  £  are  velocities  at  t  =  0.  Assuming  these 
o  o  o 


■ 
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velocities  to  be  as  i'  =  - 


9<j) 


o 


o 


9x 


we  have 


1 


<j)  =  c  /  s  dt  +  <j> 


o 


differentiating  we  have: 


4>  =  c  s 


(b .  18) 


In  one  dimension  it  has  been  shown  that  the  equation  of 
continuity  (b . 5 )  can  be  written  as: 


ii  +  i£  .  o 

9t  dX 


(b . 19) 


A  simple  extension  to  three  dimensions  results  in: 


9JL  +  ii+  o 

3t  9x  9y  9z 


(b . 20) 


Substitution  of  equations  (b.15),  (b.16),  (b . 17 )  and  (b.18) 


i 

along  with  the  approximation  <j>  =  ~  (reference  [16])  yields: 


4  -  cV*  =  0 

9t 


B.5  The  validity  of  —  >>  ^^2 

J  9x  9x 


In  order  to  check  the  validity  of  the  approximation  that 

•  9 

>>  an  examination  of  a  plane  wave  is  carried  out. 

9x  9x 

For  a  plane  wave  travelling  in  the  x  direction: 


2  2 
9  <J>  2  9  <J> 

-  —  —  r*  -  — 


9t‘ 


9x 


0 


axf3-  ... 
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The  solution  is: 

,  .  i(o)t  -  kx) 

<j)  =  A  e 

Hence 

€  =  -  =  ikA  ei(u,t  '  kx>  (b .  21) 

dX 

C  =  dt  =  el(wt  "  kx)  (b .  22) 

A  relationship  between  E,  (displacement)  and  s  (condensation) 
is  derived. 


8  £ 

6v  =  (dx  +  ^  dx) A  -  A  dx 
9x 
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v 

o 

A 

S 

A 

s 

If  the  wave  exists  in  water  at  a  depth  of  less  than  one  meter 
it  would  be  impossible  to  have  partial  pressures  as  high  as  one 
atmosphere  because  of  cavitation. 

6P  <  1.013  x  105  NT./M2 

6P  =  k's  (b . 8) 

Kinsler  and  Frey  [17]  define  a  quantity  called  the  isothermal 
bulk  modulus  (B^)  of  a  fluid  by: 


From  (b.ll)  we  have: 

2  _  kj_  „  k^_ 

=  p  =  Po 
so  that  k’  =  yB^ 

Values  for  y  and  B  are  listed  in  reference  [18]  as: 
y  =  1.004 

B  =  2.18  x  109  NT./M2 
t 


=  A  dx 


9x 


ik 


03 


A  e 


i(u3t  -  kx) 


=  ik  5 


(b.23) 


■ 
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Hen  ce  k'  =  2.19  x  109  NT. /M2  for  water. 

Now  an  upper  bound  may  be  placed  on  the  magnitude  of  s 


s|  <  -1--013  x  10-  =  .463  X  1<T4 

2.19  x  10 


S  = 

c  =  co/k 

Hence  by  equation  (b.23)  we  have 


5  =  cs 


(b.24) 


For  P  =  1  ATM.  we  have: 


ii 

9x 


=  -icks  =  (-ick)(.463  x  10  )e 


-4.  i(oot  -  kx) 


(b . 25) 


9(gs)  _  8(c§z) 


9x 


9x 


=  (-ick)(.43  x  10  )e 


-8N  2i(oot  -  kx) 


(b.26) 


g  r  ^  \ 

The  magnitude  of  t—  is  much  greater  than  that  of  v  ■■ 

OX  dX 

0*  •  v  • 

""  "t”  £  s  )  9  £ 

so  cue  approximation  - - H  t—  is  valid.  Note  also  that 

dX  dX 

•  •  ^ 

sE,  <<  £  and  —  <<  1,  two  other  approximations  used  in  this  appendix, 


B.6  Velocity  potential  in  the  horn 

An  approximate  equation  describing  the  velocity  potential 
within  the  horn  is  now  discussed.  The  approximations  used  are 


i  laq  <‘Oa  \  Si 


■ 


' 
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discussed  in  reference  [19]. 

The  equation  for  velocity  potential  which  is  used  to  describe 
wave  motion  in  the  horn  is  derived  by  assuming  that  the  waves  in  the 
horn  are  plane  and  propagate  in  the  direction  defined  by  the  axis  of 

I 

the  horn.  A  short  explanation  of  this  assumption  is  given. 

Consider  a  curvilinear  coordinate  system  as  shown  below. 


Lines  of  constant  y  are  perpendicular  to  lines  of  constant  6. 
The  distance  between  two  lines  of  constant  u  is  not  a  constant  but 
is  dependent  on  0 . 

If  the  horn  did  not  flare  rapidly  then  one  would  expect  that 
waves  would  propagate  in  the  direction  defined  by  lines  of 
constant  6.  In  this  case  deriving  an  expression  for  the  velocity 
potential  would  involve  solving  the  equation  of  continuity  and  the 


. 
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the  force  equation  for  the  element  shown  below. 


Line  of  constant  0 


S 


velocities  are  assumed  perpendicular 
to  the  cross-sectional  area 


S  =  cross-sectional  area 


For  horns  which  do  not  flare  rapidly  the  term  dp  would  be 
essentially  a  constant  for  the  element  shown  above  (reference  [20]). 

Writing  the  equations  for  continuity  and  force  in  this 
curvilinear  system  is  difficult  because  of  the  dependence  of  dp  on  0 
and  because  there  is  no  simple  expression  for  cross-sectional  area. 
To  circumvent  this  problem  a  rectangular  coordinate  system  is  used. 
The  waves  are  assumed  to  travel  in  the  x  direction  only. 


volume  element 
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For  horns  which  do  not  flare  rapidly  dx  =  dp ,  S  =  S ’  and 
9S  8S'  , 

gx  dx  -  -g  ■  dp.  Under  these  conditions  the  equations  for  continuity 
and  force  would  appear  mathematically  identical.  The  case  of  the 
plane  wave  approximation  is  relatively  easy  to  solve  and  accurately 
describes  the  amplitude  dependence  of  the  pressure  and  displacement 
velocity.  However,  the  solution  does  not  accurately  describe  the 
direction  of  propagation  or  phase  dependence  of  waves  within  the 
horn. 

In  reference  [2  ]  P.M.  Morse  claims  that  plane  wave 

r 

I  S  I 

approximation  is  valid  for  horns  satisfying  the  condition  —  <<  1. 

2m 

r 

For  the  horn  used  is  this  thesis  —  =  .03.  The  term  r  is  the 

1  2m '  s 

radius  of  the  small  end  of  the  horn. 

In  defence  of  the  approximations  used  it  can  be  said  that 
the  wave  entering  the  large  end  of  the  horn  is  plane  and  at  the 
smaller  end  the  velocity  should  be  in  the  direction  of  the  x  axis 
since  at  the  boundary  velocities  must  be  tangential  to  the  horn 
surface  which  is  almost  parallel  to  the  horn  axis.  On  the  axis  the 
velocities  are  in  the  x  direction  by  symmetry  arguments.  Recall 
too  that  the  horn  aperture  is  only  a  fraction  of  a  wavelength  in 
radius  so  that  a  velocity  in  the  y  direction  is  not  likely  to 
exist  in  the  section  of  the  horn  near  the  small  aperture. 

The  approach  used  in  this  thesis  to  solve  for  the  velocity 
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potential  of  the  horn  does  not  lead  to  an  exact  solution  but  as 
exact  solution  is,  the  author  feels,  beyond  the  scope  of  this 
thesis . 

The  equation  of  continuity  for  the  plane  wave  approximation 
is  now  derived 


£ 


S 


5  + 


3x 


dx 


S  =  cross-sectional  area 
(mx/2) 

=  S,  e 


For  details  of  horn  geometry  consult 
figure  2. 


6t 


h  p| 

3t 

Sdx  =  6t 

pis  - 

PS  +  i£^dx 

dX 


c,  ,  8S  , 

S  +  —  dx 

dX 


|f-  s  dx  =  pi  y-  dx  +  S  dx  + 

dt  dX  dX 


Ignoring  second  order  quantities  we  have: 


3pS  =  pg3S  S3(pQ 

3t  3x  3x 


Substitute  p  =  p  (1  +  s)  and  use  the  approximations 

o 

§£  <<  k  and  (£  +  si)  =  |^-  to  obtain: 

3x  3x 


|5i  +  |-  (si)  =  0 

3t  3x 


(b.27) 
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The  approximations  used  were  previously  examined  and  justified 
for  the  case  of  plane  waves.  Another  examination  will  not  be  carried 
out  for  waves  in  the  horn  but  rather  the  previous  results  will  be 
used.  Recall  equations  b.18  and  b.15. 


$  =  c  s 


These  are  substituted  in  b.27  to  obtain: 


2  8<j>  3 In  S  .  2  32<}> 

4>  =  c  —  — r -  +  c  — 

3x  3x  „  2 

3x 

2 

For  periodic  oscillations  <j>  =  -m  <f> 


+ 


_  3x 


2 


<j>  =  0 


where  m  =  flare  constant 
k  =  wave  number 


This  has  the  solution 


(b .  28) 


■ 


. 
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4>  is  the  velocity  potential  and  is  the  equation  used  in 
chapter  3.  The  velocity  components  are  found  from  4>  by 


9x 


n 

c 


9z 


These  two  terms  are  0  for  the  horn. 


The  excess  pressure  is  found  by: 
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APPENDIX  C 

ELECTRICAL  CIRCUITS 


C.l  Low  noise  preamplifier 


& 


+ 


0- 

IN 


2.2k  ^ 


2n  3906 


H  2N  3819 


1  m  S  •  22  K 


— 0 

OUT 


o 


Gated  Amplifier 
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COMPONENT  LIST  FOR 
GATED  AMPLIFIER 


Ql  = 

2N3904 

Rl 

— 

12  K 

Q2  = 

2N3904 

R2 

= 

56 

Q3  = 

2N3904 

R3 

= 

56 

Q4  = 

2N3904 

R4 

= 

68 

Q5  = 

2N3906 

R7 

= 

270 

Q6  = 

2N3904 

R8 

= 

10  K 

Q7  = 

2N3904 

R9 

= 

2.2  K 

Q8  = 

2N3904 

RIO 

= 

2.2  K 

Q9  = 

2N3904 

Rll 

= 

220 

Q10  = 

2N3904 

R13 

= 

11 

Qll  = 

2N3906 

R14 

= 

6.8  K 

Q12  = 

2N3906 

R15 

= 

330 

Q  13  = 

2N3819 

R16 

= 

3.3  K 

Cl  = 

100  pf 

R17 

= 

330 

C2  = 

250  pf 

R18 

= 

15  K 

Dl  = 

BAX13 

R19 

= 

11 

D2  = 

BAX13 

R20 

= 

52  4W 

D3  = 

BAX14 

R21 

= 

3.3  K 

D4  = 

1N4454 

R23 

= 

3.9  K 

D5  = 

1N4454 

R24 

= 

270 

R25 

22 

270 

. 

Transducer  drive  amplifier 
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! 


I 


O 


R16 
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COMPONENT  LIST  FOR 
TRANSDUCER  DRIVE  AMP 


Q1 

2 

MPSH05 

Rl 

1  K  1W 

Q2 

= 

MPSH05 

R2 

= 

470 

Q3 

= 

MPSA56 

R3 

= 

68 

Q4 

- 

MPSA06 

R4 

= 

22  K 

Q3 

= 

MPSA06 

R5 

= 

1  K 

Q6 

= 

MPSA56 

R6 

= 

330 

Q7 

= 

MPSA56 

R7 

= 

11 

Q8 

= 

MPSA06 

R8 

= 

2.2  K 

Cl 

= 

.0015  uf . 

R9 

= 

56 

C2 

= 

.005 

RIO 

= 

2.2  K 

C3 

= 

.01  uf 

Rll 

= 

68 

C4 

= 

.01  uf 

R12 

= 

68 

C5 

= 

.003  uf 

R13 

= 

120 

C6 

= 

5.aif 

R14 

= 

1.8  K  2W 

C7 

= 

.005  uf 

R15 

= 

22  K 

D1 

= 

1N4454 

R16 

= 

390 

D2 

= 

1N4454 

R17 

= 

100 

D3 

= 

1N4454 

R18 

= 

330 

D4 

— 

1N4454 

R19 

= 

11 

Amplifier  and  level  detector  circuit 
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COMPONENT  LIST  FOR  AMPLIFIER 
AND  LEVEL  DETECTOR  CIRCUIT 


Q1  = 

2N  3704 

R1  =  IK 

Q2  = 

2N  3704 

R2  =  560 

03  = 

2N  3906 

R3  =  IK 

Q4  = 

2N  3819 

R4  =  180 

Q5  = 

2N  3819 

R5  =  12K 

Q6  = 

2N  3819 

R6  =  10K 

Q7  = 

2N  3702 

R7  =  2.2K 

Q8  = 

2N  3702 

R8  =  820 

09  = 

3704 

R9  =  IK 

Q10= 

2N  3704 

R10=  1.2K 

Qll= 

2N  3704 

Rll=  2.2M 

Cl  = 

.005  uf. 

R]  2=  100K 

C2  = 

.056  uf. 

R13=  J2K 

C3  = 

360  pf. 

R14=  270 

Dl  = 

BAX  13 

R15=  2.2K 

PI  = 

IK  pot 

R16=  12K 

R17=  12K 

FIGURE  33 

Gated  amplifier  and  booster  amplifier 


FIGURE  3  4 


Amplifier  and  level  detector 
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FIGURE  35 


Electronics  rack 
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APPENDIX  D 

PROGRAM  TO  DETERMINE  DIRECTIVITY  FUNCTION 


V  R  J1PK  DTA  ;  PW ;  PC  1 ;  PC  ;  A  ;  X 1  ;  X2  ;  F 

[1]  RW+R Cl] 

[2]  RC+Rl 2] 

[3]  P+-2  4J? 

[4]  2772M«-/)27![  1  ;  ] 

[5]  AMPN+DTA12;'] 

[6]  /W  1  o  (  (  TUT  A  )x(olil30)  )  ) 

[  7  ]  /?[/•«-(  2  5 . 4  xp [ 2 ] x  0 .  5  )  *1.432 

[8]  KAW<-(  02  )  x p  /  / 

[  9]  RC+(  25. 4x/?[  2]x0  .  5  )  -M  .  48  2 

[10]  KAC+(o2)xRC 

[11]  Xl+(KAW)x(A) 

[12]  X2+(KAC)*(A) 

[13]  PW+PC+PT+\I+ 0 

[14]  G:J«-J+1 

15]  1  PltfUPPYCJim  ) 

[16]  PW+PV,  (  2xJ/f  ) 

[17]  1  Blflra?*(.Y2[I]  ) 

[18]  PC+PCA  2xJ11) 

[19]  -+G*\  (pTH7A)*I 

[20]  I+-  0 

[21]  GA-.I+I  +  l 

[22]  PP1«-P[I]  xPC 

[23]  PCl«-(  PC1+PU)  :-(  [  /(PPl+P.V)  ) 

[24]  PT+-PT ,  P<71 

[25]  -KM  x  i  (  p  P )  *  J 

[26]  V+TRTA  ,  /1/UPP  ,  P2* 

[27]  P!T/1^Z/lP^P/l[/^PPl<P71^P,7^P^Yl^X2^i  0 

[28]  »P  IS  PR ACT TOP  ATTRIBUTED  TO  COPPER' 

[29]  r >P 

[30]  ’  ' 

[31]  'HATER  APPEPATURE  RADIUS  IN  LAPP.  DA  '  \RW 

[32]  ’  » 

[33]  'COPPER  APPEPATURE  IN  LAMBDA  '  ;RX 

[34]  ’ 

ANCLE  MEASURED  AS  PEP  F' 

[35]  VH  ( ( ( p  V) - ( pTUTA ) ) , ( p  TNT A  ) ) p  V ) 

[36]  $7 


\7 


. 

. 


. 
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The  program  listed  below  evaluates  the  function: 


J-l(  x  ) 

X 


This  program  is  used  within  the  programs  listed  in 
appendices  D,  E  and  F. 


V  /•/  F3QVFXX  X 

C 1 J  X+ ( 0 , i 2  0 ) 

[  2  ]  JM+  +  /  (  (  (  ~1  )  *K  )  x  (  (  X  )  *  (  (  2  xtf )  )  )  ^  K )  x  (  »  (  K+M )  )  x  ( 

2 *(/•/+(  2xK)  )  )  )  ) 

V 


' 
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APPENDIX  E 


PROGRAM  TO  EVALUATE  INTERFERENCE  PATTERN 


V  A  PHD  I 
L1J  Al+All] 

L  2  ]  #*-26.64 

L  3  J  PH-*-AD+\  0 

L  4  J  A2+A12] 

L  5  J  .533 

L  6  J  /*-!*(  01)4180 

L'/J  0*--(  (40x(ol  )*180  )+J) 

L  8  J  START :  0*-0+J 

L  9  J  Ll*-(  (  (#x( 2oO) )*2 )+(  (  (#x(  loO)  )-A  )*2  )  )*0 . 5 
L  1 0  J  L  2*-(  (.  (#x(.  200)  )*2  )  +  (  (  loo  )  )+A  )*  2  )  )*u.  5 

L  1 1  J  (  LI -L  2  m  0(2  *0.148  )  )) 

L  12  J  ylP«-(  {  (>?1  +  M2x(  20  TH )  ))*  2  )  +  ((A2x{  lo  TH  )  )*2  )  )*u  .  b 
L  1  3  J  1  PlO/PP^UOl  ;x(  1O0)  ) 

L  1  4  J  /1PP*-J7f 
L 1 5  J  AP+APPxAP 
L  1 6  J  AD*- AD  %AP 
L 1 7  J  PH+PH , 0 

L 1 8  J  +START* iO<( 40x(oi ;*180; 

L  1  9  J  PH+-PH xl80*(Ol  ) 

L  2 0  J  /1P*->4Pt  (  \  / AD  ) 

L 2 1  J  OVT+K  (2,(p AD))p(PH,AD)  ) 

L  2  2  J  UUT+QOUT 


V 
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APPENDIX  F 

PROGRAM  WRITTEN  TO  EVALUATE  FREQUENCY  VERSUS  AMPLITUDE 


Cl] 
r  2  ] 
C3] 

[4] 

[5] 


V  R  CU  LE \AA\B \BB  ;C  ;CA;CC  ;COS \E \ETA\F  ;FO  \F1\K\L\LE  ;/'/;// ; 


PSI iQiRiRO iSIN iWiXiY 
{■!■*•  0 . 0  4  B  x  1 


F*-F*-  \  0 
LE+LE* 2. 54 


£*-148100 
F<9  *-1023000 


[  G]  RO+-0 . 8  5  */?*-/?  x  2  .  5  4 

[7]  L+C*FO 

[8]  K<- ( o 2  ):L 

[9]  //*-  r  (  (  (  (  (KxK)-(MxMi  4)  )*0. 5)xLE)*oi  ) 
[10J  LE+IlxL*  2 

[11]  £*-2  8 

[12  ]  F//*-l  0  2  0  0  0  0 

[13]  L+Q*(o'2xFN) 

[14]  CA<-H  (  Lx(  ( 02 xFiV  )*2  )  ) 

[15]  FP*-9  0  0  0  0  0 

[16]  5!T/l/?in:FP**FF+(  2  000  ) 

[17]  K+(W+o2*FR)*C 

[18]  F*-(  (FxA')-(  •‘/x.’lf-:  4  )  )*0.  5 

[19]  COS+BiK 

[20]  £\Z7/*-//*(  2  x A'  ) 

[21]  A/l**lo( F xLF  ) 

[  22  ]  BB*r(C0S*  (  20  (B*LE  )  )  )  -  (  FI//*  ( lo(  B*LE  )  )  ) 

[23]  A'*-  (  ( A'*F<9  )  *  2  )  *  (  (  (  ( Ax/?0  )  )  *  2  )  + 1  ) 

[24]  Y+KxRO  HUKxlW)*  2  )  + 1 ) 

[25]  CO/l/l+dxFF ) 

[  26  ]  P5J*-1  v  (  (  1  +  (  (  ( J/xL  )  -  (  1*  ( Jvrx0’/1  )  )  )  *  2  )  )  * 

0.5) 

[27]  ETA+COSxPSIH  (  {CCxCC)  +  iBB*X*X*BB)  )* 
0.5) 


[28]  z>r„z;zvi 

[  2  9  ]  F*-F ,  F/f 

[30]  -+S2ART*  iF/? <1100000 

[31]  UF*-((2,pF)pF,F) 

[32]  KF*-Q!'F 


V 


•  , 
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APPENDIX  H 


RADIATION  FROM  A  SIMPLE  SOURCE 


The  velocity  potential  for  a  spherical  wave  is  given  as : 


a  =  A  ei(wt  -  kr) 

r 


Hence  the  displacement  velocity  u  is  given  by: 


3d)  d) 

u  -  -  t  =  +  lk(j) 

dr  r 


p,  the  partial  pressure  is  given  by 


po8<(, 


p  =  —  =  imp  * 


Z  =  ^  =  impedance 


Z  = 


icop  <j) 
o 


itorp 


|i+ikU  1  +  ikr 


Now  consider  a  small  pulsating  sphere  of  mean  radius  a.  Let 
the  velocity  at  any  point  on  the  surface  of  the  sphere  be  given  by 


u  =  U  e 
s 


icot 


The  medium  is  in  contact  with  the  sphere  so  that  at  the 


boundary  we  have: 
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U  e 


io)t 


=  £ 


Z 


ika  +  1 

ip  ooa 
o 

If  the  radius  a  of  the  sphere  is  much  smaller  than  a  wave¬ 
length  then  ka  <<  1  so  that: 


U  e 


I0)t 


ip  toA 
o 


i(cot  -  ka) 


A  = 


2  2 
1  +  k  a 


Ua 


2 


Hence : 


ip  coa  U 
o 

P  =  - “ -  e 


P  = 


ip  cka  U 
o 


i(oot  -  kr) 
ei(wt  -  kr) 


Often  the  partial  pressure  is  written  in  the  form; 


ip  ck  .  ,  .  , 

o  ^  i(wt  -  kr) 

p  =  — : -  Q  e 

47rr  s 


where  Q  =  47ra  U 
s 


Q  is  called  the  source  strength  and  is  the  maximum  volume 
s 

flow  rate  of  fluid  through  a  surface  enclosing  the  source. 

Radiators  which  are  much  smaller  than  a  wavelength  produce 
radiation  which  is  spherical.  Suppose  that  the  radiator  is  a  small 


flat  disc  of  area  dS.  If  this  disc  is  mounted  in  a  baffle  so  that 


It 


' 
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radiation  is  confined  to  one  side  of  the  baffle  the  source 
strength  of  this  radiator  is  given  as: 


Q  =  2  U.dS 
s  —  — 

where  is  the  inner  product. 

The  volume  flow  rate  from  this  source  is  U_.  dS^  The 
factor  2  appears  because  the  radiation  is  confined  to  one  side  of 
the  baffle.  The  partial  pressure  for  this  radiator  of  area  dS^ 
which  is  much  smaller  than  a  wavelength  is  written  as: 


dp  = 


ip  ck 

o 

2iTr 


e1(u)t  '  kr)  U.dS 


■ 
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